■  Defence  Research  and  Recherche  et  developpement 

■  "  ■  Development  Canada  pour  la  defense  Canada 


DRDC I RDDC 


Analysis  of  matched  filter  mismatch  for  focusing  moving  targets  in 
multi-channel  synthetic  aperture  radar 

David  Kirkland 

DRDC  -  Ottawa  Research  Centre 


Defence  Research  and  Development  Canada 


Scientific  Report 

DRDC-RDDC-201 4-R51 

September  2014 


Canada 


Analysis  of  matched  filter  mismatch  for 
focusing  moving  targets  in  multi-channel 
synthetic  aperture  radar 

David  Kirkland 

DRDC  -  Ottawa  Research  Centre 


Defence  Research  and  Development  Canada 

Scientific  Report 
DRDC-RDDC-2014-R51 
September  2014 


©  Her  Majesty  the  Queen  in  Right  of  Canada,  as  represented  by  the  Minister  of  National 
Defence,  2014 

©  Sa  Majeste  la  Reine  (en  droit  du  Canada),  telle  que  representee  par  le  ministre  de  la 
Defense  nationale,  2014 


Abstract 


This  report  derives  the  matched  filter  necessary  for  focusing  moving  targets  in  multi¬ 
channel  Synthetic  Aperture  Radar  (SAR)  systems  after  clutter  cancellation  has  been 
performed.  The  Displaced  Phase  Centre  Antenna  (DPCA)  is  utilized  to  achieve 
the  clutter  cancellation  in  a  two  channel  SAR  system.  After  deriving  the  matched 
filter,  the  tolerance  of  the  filter  is  then  analyzed  for  mismatches  against  errors  in  the 
moving  target’s  position  and  velocity  components.  The  tolerances  are  quantified  for 
two  exemplar  SAR  systems;  an  X-band  airborne  SAR  and  a  C-band  satellite  SAR. 
The  analysis  also  reveals  when  simplified  versions  of  the  matched  filter  can  be  used. 

Significance  for  defence  and  security 


The  ability  to  generate  a  focused  image  of  a  moving  target  buried  in  stationary  clutter 
can  provide  increased  intelligence  and  can  improve  the  performance  of  automatic 
target  recognition  (ATR)  algorithms.  Moving  Target  Indication  (MTI)  techniques  are 
often  focused  at  detecting  slow  moving  targets.  In  SAR  images  these  slow  moving 
targets  are  fairly  well-focused,  although  they  may  be  displaced  in  the  image  from 
their  true  location.  The  use  of  Space-Time  Adaptive  Processing  (STAP)  techniques  in 
combination  with  SAR  imaging  has  been  discussed  in  the  literature  already.  However, 
the  exact  nature  of  the  resulting  moving  target  signature  and  how  to  produce  a  focused 
image  of  the  moving  target  has  not  been  extensively  studied.  This  report  is  aimed 
addressing  this  shortcoming.  Additionally,  the  parameters  extracted  from  focusing 
the  moving  target  can  utilized  in  tracking  algorithms. 

This  report  derives  the  moving  target’s  velocity  tolerances  necessary  to  produce  a 
focused  image  of  the  moving  target.  These  velocity  tolerances  are  calculated  for  a 
typical  airborne  X-band  SAR  system  and  for  a  satellite  C-band  SAR  system  with 
parameters  similar  to  RADARS  AT- 2. 
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Resume 


Le  present  rapport  indique  les  calculs  relatifs  au  filtre  adapte  necessaire  pour  focaliser 
des  objectifs  mobiles  dans  les  systemes  radar  a  synthese  d’ouverture  (SAR)  multica- 
naux  apres  avoir  climine  les  echos  parasites  (clutter).  L’antenne  a  centre  de  phase 
deplace  (DPCA)  est  utilisee  pour  eliminer  le  clutter  dans  un  systeme  SAR  a  deux 
canaux.  Une  fois  que  les  caracteristiques  du  filtre  adapte  sont  calculees,  la  tolerance 
de  ce  filtre  est  analysee  pour  trouver  les  defauts  d’adaptation  par  rapport  aux  erreurs 
dans  les  composantes  de  position  et  de  vitesse  de  hobjectif  mobile.  La  tolerance  est 
calculee  pour  deux  systemes  SAR  de  reference  :  un  SAR  aeroporte  en  bande  X  et  un 
SAR  satellitaire  en  bande  C.  L’analyse  montre  aussi  a  quel  moment  on  peut  utiliscr 
des  versions  simplifies  du  filtre  adapte. 

Importance  pour  la  defense  et  la  securite 


La  capacite  de  produire  une  image  focalisee  d’un  objectif  mobile  cache  dans  du  clut¬ 
ter  stationnaire  peut  fournir  des  renseignements  supplementaires  et  ameliorer  les  per¬ 
formances  des  algorithmes  de  reconnaissance  automatique  des  objectifs  (ATR).  Les 
techniques  de  visualisation  des  cibles  mobiles  (VCM)  visent  souvent  a  detecter  des 
objectifs  mobiles  lents.  Dans  les  images  SAR,  ces  objectifs  mobiles  lents  sont  assez 
bien  Realises,  meme  s’ils  sont  parfois  deplaces  dans  Limage  par  rapport  a  leur  position 
reelle.  L’utilisation  conjointe  des  techniques  de  traitement  spatio-temporel  adaptatif 
(STAP)  et  de  l’imagerie  SAR  a  deja  fait  l’objet  d’etudes  .  Par  contre,  la  nature  exacte 
de  la  signature  resultante  d’un  objectif  mobile  et  la  production  d’une  image  focalisee 
de  hobjectif  mobile  n’ont  pas  ete  etudiees  a  fond.  Le  present  rapport  vise  a  corriger 
cette  lacune.  De  plus,  les  parametres  extraits  de  la  Realisation  de  hobjectif  mobile 
peuvent  etre  utilises  dans  les  algorithmes  de  poursuite. 

Le  present  rapport  indique  les  calculs  relatifs  aux  tolerances  de  vitesse  d’un  objectif 
mobile  qui  sont  necessaires  pour  produire  une  image  focalisee  de  cet  objectif.  Ces 
tolerances  de  vitesse  sont  calculees  pour  un  systeme  SAR  aeroporte  type  en  bande  X 
et  un  systeme  SAR  satellitaire  en  bande  C  ayant  des  parametres  similaires  a  ceux  de 
RADARS  AT  2. 


DRDC-RDDC-2014-R51 


Table  of  contents 

Abstract .  i 

Significance  for  defence  and  security  .  i 

Resume .  ii 

Importance  pour  la  defense  et  la  securite .  ii 

Table  of  contents .  iii 

List  of  figures .  iv 

1  Introduction .  1 

1.1  Mathematical  notation  .  1 

1.2  Signal  model .  3 

1.3  DPCA  imaging  geometry .  4 

1.4  Moving  target .  5 

2  Matched  filter  analysis .  10 

2.1  DPCA  matched  filter .  10 

2.2  DPCA  matched  filter  mismatch  .  15 

3  Summary  and  conclusions .  23 

References .  25 

Annex  A:  Taylor  series  approximations  for  cos^  and  sin  02  .  27 

A.l  Approximation  for  cos  02  (w) .  27 

A. 2  Approximation  for  sin  02  (u) .  28 


DRDC-RDDC-2014-R51 


List  of  figures 


Figure  1:  SAR  Imaging  Geometry .  3 

Figure  2:  DPCA  Geometry. .  5 


IV 


DRDC-RDDC-2014-R51 


1  Introduction 


The  effect  of  a  moving  target  in  Synthetic  Aperture  Radar  (SAR)  imagery  was  first 
analyzed  in  [1],  Subsequently,  algorithms  were  developed  to  focus  the  moving  target 
within  the  SAR  image  for  single  channel  systems  [2,  3,  4,  5].  Unfortunately,  it  is  not 
possible  to  simultaneously  focus  the  moving  target  and  the  stationary  background  in  a 
single  channel  system.  Additionally,  it  is  often  only  possible  to  focus  the  moving  target 
response  when  the  stationary  background  clutter  is  relatively  weak  in  comparison  with 
the  moving  target,  since  the  motion  parameters  necessary  for  focusing  the  moving 
target  have  to  be  estimated. 

Ground  Moving  Target  Indication  (GMTI)  systems  utilize  multiple  channels  to  effec¬ 
tively  cancel  out  the  stationary  clutter  and  detect  a  moving  target  [6].  A  key  metric 
of  interest  for  GMTI  applications  is  the  minimum  detectable  velocity  (MDV).  Faster 
moving  targets  are  usually  easier  to  detect  due  to  the  increased  Doppler  separation 
between  the  moving  target  signal  and  the  stationary  background  clutter  [7].  SAR 
processing  for  stationary  targets  will  produce  a  fairly  well-focused  SAR  image  for 
slow  moving  targets.  In  terms  of  focusing  moving  targets,  the  main  interest  lies  in 
the  higher  velocity  targets,  which  produce  smearing  and  range  migration  effects  in 
the  resulting  uncompensated  SAR  imagery  [1], 

This  report  conducts  a  theoretical  analysis  of  the  use  of  a  multi-channel  SAR  system 
in  a  Displaced  Phase  Centre  Antenna  (DPCA)  configuration  to  allow  the  focusing 
of  the  moving  target  signature  when  it  is  buried  in  the  stationary  clutter.  Section 
1.2  develops  the  signal  model  for  the  Spotlight  SAR  model.  The  Spotlight  SAR 
configuration  is  used  since  it  allows  the  formation  of  higher  resolution  imagery  in 
comparison  to  Stripmap  imagery  operating  with  the  same  parameters.  Section  1.3 
describes  the  DPCA  imaging  geometry.  Section  1.4  describes  the  signal  models  for  the 
stationary  clutter,  the  moving  target  signature  and  the  resulting  signal  after  DPCA 
clutter  suppression.  Section  2  examines  the  output  of  the  matched  filter  and  the 
sensitivities  to  the  mismatch  of  the  moving  target  motion  parameters. 

1.1  Mathematical  notation 

In  this  report  we  deal  with  signals  in  both  temporal  and  spatial  dimensions.  For 
a  temporal  signal,  the  variables  t  and  c o  denote  time  and  its  Fourier  counterpart 
temporal  frequency  respectively.  The  units  of  t  and  c o  are  seconds  and  radians/second 
respectively.  The  temporal  Fourier  transform  of  a  signal  s(t)  is  denoted  as  S(u), 
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where: 


S(u)  =  Ft[s(t)] 

oo 

=  j  s(t)  exp(-juit)  dt. 

—  OO 

The  inverse  temporal  Fourier  transform  of  a  signal  S(oo)  is  identified  as: 

=  K'  PH] 

1  °° 

—  —  S(oj)  exp(jujt)  dco. 

Z7T  J 

— OO 

For  convenience,  we  define  the  wavenumber,  k,  as: 

k=-, 

c 

where  c  is  the  speed  of  light.  The  unit  of  k  is  radians/meter. 


(1) 


(2) 


(3) 


For  signals  in  the  spatial  domain  u,  the  Fourier  counterpart  is  spatial  frequency  ku, 
which  is  also  called  the  wavenumber.  The  units  of  u  and  ku  are  meters  and  radians 
per  meter  respectively.  The  forward  Fourier  transform  of  a  spatial  domain  signal  g{u) 
is  denoted  as  G(ku),  where: 


G(ku)  =  Tu[g{u)} 

OO 

=  j  g(u)  exp(-jkuu)  du. 

—  OO 

The  inverse  spatial  Fourier  transform  of  a  signal  G(ku )  is  identified  as: 

a(u)  =  ny  [g(m 

=  —  J  G(ku)exp(jkuu)  dku. 


(4) 


(5) 


For  convenience,  we  will  use  the  symbol  to  denote  a  Fourier  transform  pair,  that 
is: 


s(t)  S(uj), 

(6) 

g (u)  #  G(ku). 

(7) 

2 


DRDC-RDDC-2014-R51 


1.2  Signal  model 


The  data  collection  geometry  for  the  radar  is  shown  in  Figure  1.  The  instantaneous 
radar  position  is  given  by  (0,  m),  where  u  denotes  the  position  of  the  radar  platform 
in  the  cross-range  direction.  Xc  and  Yc  denote  the  centre  of  the  spotlight  scene  in  the 
range  and  cross-range  directions  respectively.  The  instantaneous  position  of  a  moving 
target  is  given  by  (Xc  +  x0  +  vxu,  Yc  +  y0  +  vyu),  where  x0  and  y0  denote  the  initial 
offset  position  of  the  moving  target  from  the  scene  centre  in  the  range  and  cross-range 
directions  respectively.  vx  and  vy  denote  the  moving  target’s  velocities  in  the  range 
and  cross-range  directions  respectively.  For  convenience  the  velocity  components  of 
the  moving  target  have  been  normalized  by  the  radar  platform  velocity. 


(Xc  +  x0  +  vxu,  Yc  +  y0  +  vyu) 


Figure  1:  SAR  Imaging  Geometry. 


After  range  compression,  basebanding  and  motion  compensation  to  a  point,  the  re¬ 
ceived  radar  signal  for  the  moving  target,  S(oj,u),  in  the  frequency /cross-range  do¬ 
main  can  be  written  as: 


S(uj,u)  =  P(uj)exp(j(j)(uj,u)) ,  (8) 

where: 

4>(lu,  u )  =  —2 (k  +  k0)RA(u ),  (9) 

Ra(u)  =  R(u)  —  Rref(u ),  (10) 

R(u)  =  \J (Xc  +  x0  +  vxu )2  +  (Yc  +  y0  +  vyu  -  u )2,  (11) 
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Rrefiu)  =  yj(. xc)2  +  (Yc  -  u )2,  (12) 

where  P{oj)  denotes  the  Fourier  transform  of  the  pulse  compressed  envelope,  k0  is  the 
wavenumber  corresponding  to  the  radar  centre  frequency  /0,  c  is  the  speed  of  light, 
k  —  -  is  the  fast-time  wavenumber  (relative  to  k0),  and  cj)(uj,u )  denotes  the  phase 
response  of  the  moving  target.  Rref{u )  is  the  reference  range,  and  corresponds  to  the 
instantaneous  distance  between  the  radar  platform  and  the  imaging  scene  centre.  In 
the  subsequent  section  dealing  with  DPCA,  each  channel  will  have  a  response  similar 
to  that  shown  in  (8)  but  each  channel  will  have  an  unique  instantaneous  range  R(u) 
and  reference  range  Rref(u). 

This  development  of  the  signal  model  has  neglected  to  include  the  antenna  beam 
pattern.  The  impact  of  the  antenna  beam  pattern  on  the  signal  model  developed 
in  (8),  is  to  introduce  an  amplitude  modulation  in  the  cross-range  direction.  For 
the  purpose  of  image  formation  we  are  mostly  interested  in  the  phase  portion  of  the 
response  and  therefore  the  antenna  pattern  has  been  omitted  from  (8)  and  from  the 
subsequent  analysis.  In  order  for  a  DPCA  system  to  achieve  clutter  cancellation 
the  two-way  beam  patterns  for  the  Transmit /Receive  pairs  must  be  identical.  When 
the  beam  patterns  are  not  identical  then  channel  balancing  procedure  is  required  to 
achieve  effective  cancellation. 

1.3  DPCA  imaging  geometry 

The  Displaced  Phase  Centre  Array  (DPCA)  geometry  typically  involves  a  single  trans¬ 
mitter  and  two  receivers.  The  DPCA  condition  implies  that  the  array  moves  a  specific 
distance  between  consecutive  pulses  such  that  adjacent  pulses  at  the  two  receivers 
can  be  subtracted  to  remove  the  stationary  clutter. 

In  this  report  a  modified  array  geometry  is  utilized  to  ease  the  burden  of  the  mathe¬ 
matical  analysis  that  follows  which  examines  the  sensitivity  to  the  mismatch  of  focus¬ 
ing  parameters.  In  this  geometry  we  assume  that  there  are  two  transmitter /receiver 
pairs  denoted  as  Tx\ ,  Rx i,  Tx 2,  Rx 2.  Tx\  and  Rx\  are  collocated  and  Tx 2  and 
Rx 2  are  collocated.  Tx±  and  Tx 2  are  separated  by  a  distance  d.  The  SAR  system 
is  assumed  to  move  the  distance  d  between  consecutive  pulses  i.e. ,  meets  the  DPCA 
condition.  This  geometry  is  illustrated  in  Figure  2. 

The  position  of  Tx\/Rx\  is  given  by  (0,  u),  and  the  position  of  Tx2/Rx 2  is  given  by 
(0,  u  —  d).  The  instantaneous  range  between  Tx\jRx\  and  the  stationary  target  point 
(Xc  +  Xo,  Yc  +  y0 )  is  denoted  as  R\[u)  and  is  given  by: 

Ri(u)  =  sj {Xc  +  x0)  +  (Yc  +  yo~  u )2.  (13) 
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Similarly,  the  instantaneous  range  between  Tx2/ Rx2  and  target  point  (Xc+Xo,  Yc+y0) 
is  denoted  as  R2(u)  and  is  given  by: 

R2{u)  =  \] (Xc  +  x0 )2  +  (Yc  +  y0  -  u  +  d )2.  (14) 

It  is  straightforward  to  see  that  R\(u  —  d)  =  R2{u ),  and  therefore  the  DPCA  condition 
is  met  when  the  inter-pulse  distance  is  equal  to  d. 

1.4  Moving  target 

We  now  develop  the  range  equations  for  a  moving  target.  Consider  a  target  with 
velocities  vx  and  vy  in  the  range  and  cross-range  directions  respectively.  As  before, 
the  velocities  have  been  normalized  by  the  radar  platform  velocity.  If  the  initial 
position  of  the  target  is  (Xc  +  £0,  Yc  +  yo),  then  instantaneous  position  of  the  moving 
target  is  given  by  (Xc  +  x0  +  vxu,  Yc  +  y0  +  vyu )  and  the  instantaneous  range  between 
Txi/Rxi  and  moving  target  is  given  by: 

Ri(u)  =  yj (Xc  +  x0  +  vxu )2  +  (Yc  +  yo  +  vyu  -  u )2.  (15) 

Similarly,  the  instantaneous  range  between  Tx2/Rx2  and  moving  target  is  given  by: 

R2(u)  =  sj  (Xc  +  x0  +  vxu )2  +  (Yc  +  y0  +  vyu-u  +  d )2.  (16) 

For  convenience  we  dehne  two  reference  ranges.  These  reference  ranges  are  used  in  the 
motion  compensation  to  a  point  procedure  necessary  for  Spotlight  image  formation. 
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The  reference  range  Rrefi(u)  denotes  the  distance  between  the  Tx±  and  the  scene 
centre  (Xc,  Yc),  and  is  given  by: 

RrefM  =  ^(xcy  +  {Yc-uy.  (17) 

The  second  reference  range,  Rref2(u)  denotes  the  distance  between  the  Tx 2  and  the 
scene  centre  (XC,YC),  and  is  given  by: 

Rref2(u)  =  s/(Xcy  +  {Yc-U  +  d) 2.  (18) 

Again  it  is  straightforward  to  verify  that  Rref  \  (u  —  d)  =  Rref2(u). 

In  general  the  received  signal  after  range  compression  in  the  fast-time  frequency  / 
cross-range  domain  is  denoted  by  Srx(k,u),  where: 

Srx(k,  u )  =  P(u)  exp(—j2(k  +  k0)R(u)),  (19) 

where  k  is  the  fast-time  wavenumber,  k0  =  2irf0/c  is  the  wavenumber  corresponding 
to  the  radar  centre  frequency  /o,  P{oj)  is  the  pulse  compression  response.  After 
performing  the  motion  compensation  to  a  point,  the  resulting  signal,  S(k,u)  is  given 
by: 


S(k,  u )  =  P{oj)  exp(— j2(k  +  ko)(R(u)  —  Rref{u ))),  (20) 


where  Rref(u )  is  the  instantaneous  distance  from  the  platform  to  the  scene  centre. 

The  signal  Si(k,u)  from  the  Tx\jRx\  pair  is  then  given  by: 

Si(k,u)  =  P(uj)exp(-j2(k  +  k0)(Ri(u)  -  Rrefl(u))) 

=  P(u)  exp(—  j2(k  +  ko)ARi(u)), 


where  ARi(u)  =  Ri(u)  —  Rrefi(u).  Similarly,  the  signal  S2(k,u )  from  the  Tx2/Rx2 
pair  is  then  given  by: 

S2(k,u )  =  P(u)exp(-j2(k  +  k0)(R2(u)  -  Rref2(u))) 

—  P(co)  exp(—j2(k  +  k0)AR2(u)), 


where  A R2{u)  =  R2(u)—Rref2(u).  It  should  be  noted  that  the  removal  of  the  reference 
ranges,  Rrefi(u)  and  Rref2(u),  in  this  manner  essentially  applies  a  deramping  function 
in  the  cross-range  dimension  of  the  data.  The  result  of  this  deramping  procedure  is 
that  the  reference  point  or  scene  center  is  focused.  This  procedure  is  often  referred 
to  as  motion  compensation  to  a  point  and  is  often  used  in  the  Polar  Format  imaging 
algorithm  for  Spotlight  SAR  imaging  [8,  9]. 
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The  result  of  the  DPCA  cancellation,  Sc(k,u)  is  given  by: 

Sc{k,u)  =  Si(A;,u-  d)  -  S2{k,u) 

=  P(lu)  [exp(—j2(k  +  k0)ARi(u  -  d))  -  exp (—j(k  +  k0)AR2(u))] . 

From  the  previous  properties  of  the  range  equations  it  can  be  seen  that  the  stationary 
targets  will  be  removed  from  the  signal  Sc(k,  u),  while  the  moving  targets  will  remain. 
Expanding  A Ri(u),  and  AR2(u)  in  a  Taylor  series  w.r.t.  the  variables  x0,  y0,  vx,  vy 
about  the  point  xq  —  yo  —  vx  —  vy  —  0  gives  (see  [4]): 


A  p  ,  x  Xc(x o  +  vxu )  (Yc  -  u)(y0  +  vyu ) 

U  i2re/l(ti)  RrefM 

(24) 

Xc(x0  +  vxu)  (Yc-u  +  d)(y0  +  vyu) 
l\R2[u)  «  +  fn.\ 

(25) 

=  cos  O2(u)(x0  +  vxu)  +  sin  92(u)(y0  +  vyu), 

where: 

cos  e2{u)  = 

Rref2  (v*) 

(26) 

and 

■  a  t  \  (Yc  ~  u  +  d) 
sm  02(u)=  . 

f^ref2  \^J 

(27) 

It  should  be  noted  that: 

Rrefliy*  d)  Rref2{^)  5 

(28) 

and 


A  Ri(u  —  d ) 


Xc(xp  +  vx(u-d))  (Yc-u  +  d)(y0  +  vy(u  -  d)) 
Pref2 (w)  Rref2(td) 

cos  62{u){xq  +  vxu  —  vxd )  +  sin  d2(u)(yo  +  vyu  —  vyd). 


(29) 


It  has  been  shown  previously  how  squint  Spotlight  mode  operation  can  be  treated  as 
a  rotated  geometry  of  the  broadside  Spotlight  mode  [4,  5,  10].  For  convenience  we 
then  examine  the  broadside  case  where  02  (0)  =  0,  without  any  loss  of  generality.  In 
Appendix  A  the  following  Taylor  series  approximations  are  developed: 

„  ,  .  d  u2 

coSe2iu)*i  +  -u-  — 
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and 


■  a  /  \  d  u 

Sin  02  w  ~  -  TT- 

•^C  ^  c 

Applying  these  approximations  to  (25)  and  (29)  gives: 


+  2X1 

d 


d  u2  \  .  /  d  u  \  .  ,, 

(®o  +  vxu  -  vxd)  +  —  -  —  J  (2/0  +  vyu  -  vyd) 


d2 


=  x0  +  —y0  -  dvx  -  —vy  +  \vx  + 


VX  XCf 

dx  o  d2vx  yo  2  dv. 


X'y'J  Xr  V  '  VX  '  A?  A?  X  '  Xr. 


u 


3dvx  —  xo  t’: 
2A2  X 


y  „,2  Ux  „,3 


u 


2xf  ' 


(30) 


and 

AR2(u ) 


(  d  u2  \  .  /  d 

'  +  xf“2xd  + 


u 


(yo  +  VyU ) 


d 


dXn 


dvv 


i  i  i  i  1/0  Vo  .  ,  . 

I"  +  x90+  B'+xF"x:+x  "  + 


vy ■  »' 

/cl  dvx-x0  Vy 


2A2 


X 


2  ^  3 

W  XT7T7  ^  , 


2A2 
(31) 


Examination  of  the  components  of  Ai2i(w  —  d),  and  AR2(u)  reveals  that  they  are 
quite  similar  in  structure,  but  the  presence  of  the  moving  target  introduces  some 
additional  terms.  In  fact,  it  is  these  terms  which  enables  DPCA  and  Space  Time 
Adaptive  Processing  (STAP)  techniques  to  detect  moving  targets  in  clutter.  There  is 
a  constant  difference  of  —vxd—  ^~vy  between  ARi(u  —  d)  and  A R2(u).  The  dominate 
component  of  this  offset  will  usually  be  the  —vxd  term,  since  d  <C  Xc.  The  linear 
components  of  AR\(u  —  d)  and  AR2(u)  differ  by  the  term  —jAu.  An  extra  term, 
given  by  appears  in  the  quadratic  component  of  A Ri(u  —  d).  The  cubic  terms 
of  ARi(u  —  d)  and  A R2(u)  are  identical. 

In  deriving  these  equations  a  flat-earth  model  has  been  used  to  simply  the  nature 
of  the  calculations.  For  an  airborne  geometry  this  assumption  is  usually  adequate. 
For  spaceborne  SAR  geometries  additional  adjustments  are  required  to  accommodate 
the  curvature  of  the  Earth  and  the  gravitational  force  the  Earth  contributes  to  the 
equations  of  motion  of  the  SAR  platform.  In  [11,  12]  it  is  shown  that  for  a  single 
channel  system  the  SAR  equations  need  only  a  slight  modification  by  defining  an 
“effective”  satellite  velocity.  The  SAR  equations  for  a  Multi-Channel  Stripmap  mode 
geometry  are  developed  in  [13].  In  the  analysis  of  this  report  the  reference  range  from 
the  satellite  antennas  to  a  stationary  point  on  the  Earth  is  removed  from  the  signal 
model  (see  (21)  and  (22)).  The  reference  ranges  are  exact  i.e. ,  no  approximation  has 
been  made  -  however,  in  the  satellite  case  the  reference  ranges  will  no  longer  be  of  the 
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form  given  in  (17)  and  (18).  The  result  of  removing  these  reference  ranges  is  the  the 
reference  point  is  focused  in  the  SAR  signal.  This  is  the  first  processing  step  in  the 
Polar  format  algorithm  for  Spotlight  mode  imaging.  Therefore,  the  resulting  Taylor 
series  which  are  developed  for  ARi  and  A R2  will  still  be  applicable  for  the  satellite 
geometry. 
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2  Matched  filter  analysis 


In  section  1.4  we  developed  the  expression  for  the  moving  target  signal  after  the 
cancellation  of  the  stationary  clutter.  To  form  a  focused  image  of  the  moving  target 
a  matched  filter  needs  to  be  applied  to  the  resulting  moving  target  signature.  Section 
2.1  develops  the  necessary  matched  filter  and  derives  the  range  of  parameters  over 
which  the  output  of  the  matched  filter  can  be  expected  to  produce  focused  targets. 
Section  2.2  analyzes  the  output  of  the  matched  filter  when  the  filter  is  mismatched 
to  the  moving  target’s  motion  parameters. 

2.1  DPCA  matched  filter 

After  DPCA  cancellation,  a  focused  image  of  the  moving  target  can  be  formed  by  ap¬ 
plying  a  matched  filter  to  (23).  In  some  sense,  this  is  similar  to  time  domain  matched 
filter  processing.  It  forms  an  image  of  desired  point  and  assumes  an  invariance  region 
around  that  point  i.e  the  centre  point  is  focused  but  points  further  away  suffer  from 
some  defocusing.  An  alternative,  is  to  use  Polar  format  processing  and  apply  the 
matched  filter  afterwards. 

Based  on  the  output  of  the  stationary  clutter  cancellation  in  (23)  the  matched  filter 
M(k,  u)  for  focusing  a  moving  target  is  given  by: 

M(k,  u)  =  S*(k,  u) 

=  exp(j2(k  +  k0)ARi(u  —  d))  —  exp(j2(/c  +  k0)AR2(u)). 

It  should  be  noted  that  the  P{oj)  factor  of  the  signal  is  not  included  in  (32)  since  it  is 
assumed  that  pulse  compression  has  already  been  performed.  The  result  of  applying 
the  matched  filter  to  (23)  is  given  by: 

Sc(k,  u)M(k ,  u)  =  P{oj)  [exp(jO)  +  exp(jO)  —  exp(—  j2{k  +  k0)(ARi(u  —  d)  —  A R2(u))) 

—  exp(—j'2(k  +  ko)(AR2(u)  —  ARi(u  —  d)))] 

=  P(u>)  [2  —  exp(— j2(k  +  ko)(ARi(u  —  d)  —  A R2(u))) 

-  exp(j2(/c  +  /c0)(A-Ri(u  -  d)  —  Ai?2(u)))] 

=  P(u>)  [2  —  2  cos(2(/c  +  ko)(ARi(u  —  d)  —  AR2(u)))] . 

(33) 
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In  broadside  mode  the  expression  for  ARi(u  —  d)  —  AR2(u)  simplifies  to: 

.  „  ,  d  ,  d2  (  dxn  d2vx  yn 

ARx(u  -  d)  -  AR2(u)  =  x0  +  —y0  -  dvx  -  —  vy  +  (  vx  +  —  -  + 

c  c  \  c  c  * 


3dna;  —  xq  Wy 


2X2 


X  I  U  2X2 


®  „,3 


U 


d  (  dx  o  2/o  ,  dvy 

x0  +  —2/o  +1^  +  —  -  —  +  —  )  tt 


Xc2  Xc  Xc 


x0  -  2dnx  n,,  \  2 
+  -rf  I  u 


2X2 

d2 


X, 


s  u3 


2X2 


dv,,  d2vr 


=  -dvx  -  —vv  +  -/  - 


xc  y  '  v^c  x2 

dvy  vxd  2 


u  + 


2X 


where  it  is  assumed  that  d  -C  Xc.  Substituting  (34)  into  (33)  gives: 


2  dvy 

XT 


u 


(34) 


Sc(k,u)M(k,u)  =  P(uj) 


2  —  2  cos  (  2{k  +  fco)  (  —  +  ~rru  + 


vxd 


X,  2X2 


^u,  2 

U 


(35) 


A  focused  image  of  the  moving  target  can  then  by  obtained  by  taking  the  two  di¬ 
mensional  Inverse  Fourier  transform  of  (35).  The  2 P(oS)  term  in  (35)  is  the  baseline 
coherent  gain  of  the  matched  filter.  The  cosine-based  term  will  introduce  spurious 
components  into  the  SAR  imagery  which  are  dependent  on  the  moving  target’s  ve¬ 
locity  parameters  and  the  distance  between  the  radar  platforms  i.e. ,  d. 

The  koVxd  terms  in  the  cosine  portion  of  (35)  will  modulate  the  quiescent  coherent 
gain  i.e.,  it  can  add  either  constructively  or  destructively  to  the  constant  coherent 
gain  provided  by  the  first  term  of  (35).  The  destructive  interference  leads  to  blind 
velocities  and  will  be  discussed  further  later  on.  It  should  be  noted  that  the  dominant 
portion  of  this  term  is  due  to  the  vx  velocity  since  d  is  usually  much  smaller  than 
Xc.  The  kvxd  terms  in  the  cosine  portion  of  (35)  will  result  in  duplicate  targets 
appearing  in  the  range  direction  in  the  resulting  SAR  imagery.  The  (k  +  k0)  ^-u  term 
in  the  cosine  portion  of  (35)  results  in  range  walk  of  the  moving  target  and  duplicate 
targets  appearing  in  the  cross-range  dimension  of  the  resulting  SAR  imagery.  The 
{k  +  ko )  u2  term  in  the  cosine  portion  of  (35)  will  produce  range  curvature  and 
smearing  in  the  cross-range  dimension  of  the  moving  target  in  the  corresponding 
SAR  imagery. 

An  alternative  viewpoint  is  to  examine  the  effects  of  the  target’s  velocity  components 
on  the  output  of  the  matched  filter  operation.  The  vy  velocity  component  of  the 
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moving  target  will  introduce  duplicate  targets  in  the  cross-range  and  introduce  range 
walk.  The  vx  velocity  component  of  the  moving  target  will  introduce  duplicate  peaks 
in  the  range  direction  (due  to  the  vxd  term),  cross-range  smear  and  a  range  curvature 
component. 


The  output  of  the  matched  filter  can  be  used  to  form  a  focused  image  of  the  moving 
target  if  the  linear  and  quadratic  components  of  the  cosine  term  in  (35)  are  limited. 
These  limits  will  determine  the  maximum  target  velocities  that  can  be  imaged  with 
the  matched  filter  under  the  operating  conditions  of  the  radar  platform  i.e. ,  aperture 
length,  antenna  separation,  and  operating  frequency.  The  limits  for  the  magnitudes 
of  the  target  velocities  can  be  established  by  limiting  the  range  migration  (range  walk 
and  range  curvature)  to  be  less  than  a  range  cell,  limiting  duplicate  targets  to  be  less 
than  a  resolution  cell  distance  away  (range  and  cross-range),  and  limiting  quadratic 
phase  errors  to  be  less  than  7t/4  [11,  8].  This  results  in  the  following  five  criteria: 


1.  To  limit  the  range  walk  the  shift  should  be  less  than  one  range  cell  at  the 
extreme  limit  of  the  aperture  i.e.,  umax : 


2  kn 


\vy  \  d 


X 


l^max  27T 

t rX 


(36) 


l”»l  < 


u 


h  d 

Tx,m.a.x  w 


max  ^ max 


2.  To  limit  the  duplicate  targets  in  the  cross-range  dimension  the  cross-range  fre¬ 
quency  component  should  be  less  than  the  cross-range  resolution  i.e.,  27t/L, 
where  L  =  umax  ~  umin.  This  results  in  the  following  constraint: 


2/cn 


vyd  |  27T 


X. 


\vJ  < 


7TX 

Lknd 


(37) 


3.  To  limit  the  cross-range  smear  we  limit  the  phase  variation  at  the  limit  of  the 
aperture  i.e.,  umax : 


2kr 


luj  d  n  7T 

- ?/  <T  — 

2X2  max  4 


7rX 


(38) 


<  kr)dui, 


vQU^max 

4.  To  limit  the  range  curvature  the  shift  should  be  less  than  one  range  cell  at  the 
extreme  limit  of  the  aperture  i.e.,  umax : 

\vT\  d 


2  kn 


2A? 


InJ  <  2n 


X2 


(39) 


k  dn 2 
rhmaxUjU'max 
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5.  The  limit  the  duplicate  targets  in  the  range  direction  we  limit  the  shift  to  be 
less  than  one  range  cell,  which  produces  the  following  constraint: 


2 krnax  \Vr\  d  <C  27T 


7 r 


U  A 


(40) 


It  should  be  noted  that  all  of  these  criterion  are  inversely  proportional  to  the  phase 
centre  separation.  Therefore  a  smaller  phase  centre  separation  increases  the  range  of 
velocities  over  which  the  matched  filter  will  focus  the  moving  target. 


A  further  restriction  occurs  if  we  consider  the  operation  of  the  system  in  conditions 
where  the  previous  five  criterion  are  all  met.  In  this  case  the  terms  dependent  on  u 
become  negligible  and  (35)  is  then  reduced  to: 

Sc(k,  u)M(k,  u )  =  P(uj)  [2  —  2  cos  (2 (k  +  ho)  (— vxd))] .  (41) 


In  this  form  it  is  straightforward  to  see  that  the  DPCA  approach  is  prone  to  blind 
speeds  when: 


v. 


x  =  Nik 

at 


(42) 


0 


where  N  is  an  integer.  At  these  velocities  the  output  of  the  canceler  is  almost  zero 
over  the  bandwidth  of  the  radar. 


To  illustrate  the  expected  maximum  velocities  for  typical  radar  systems,  we  have 
considered  an  X-band  airborne  radar  system  and  a  C-band  satellite  radar  system. 
The  parameters  for  the  airborne  and  satellite  radar  systems  are  given  in  Table  1 
and  Table  2  respectively.  Note  that  the  both  the  theoretical  cross-range  resolution 
and  pulse  bandwidth  are  kept  constant  for  both  systems.  The  pulse  bandwidth  of 
k0  =  4.1867  rad/m  corresponds  to  200  MHz. 


Table  1:  Airborne  radar  parameters. 


Parameter 

Value 

k0 

188.5  rad/m 

h 

rvmax 

4.1867  rad/m 

d 

1  m 

^ max 

511  m 

L 

1022  m 

xc 

20,000  m 

The  calculated  velocity  tolerances  for  the  airborne  and  satellite  radar  systems  are 
shown  in  Table  3  and  Table  4  respectively.  It  should  also  be  noted  that  the  satel¬ 
lite  velocity  is  on  the  order  7100  m/s  and  therefore  the  aperture  illumination  time 
corresponds  to  73796/7100  =  10.4  seconds. 
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Table  2:  Satellite  radar  parameters. 


Parameter 

Value 

k0 

111  rad/m 

h 

•^max 

4.1867  rad/m 

d 

7.5  m 

V'  max 

36898  m 

L 

73796  m 

Xc 

850  km 

Table  3:  Maximum  velocity  constraints  for  an  airborne  radar  with  the  operating 
parameters  given  in  Table  1.  For  the  blind  velocities,  N  denotes  any  integer. 
Velocities  are  normalized  by  the  platform  velocity. 


Term 

Limit 

Range  Walk 

\vy\  <  29.37 

Cross-Range  Duplicates 

\vy\  <  0.3262 

Range  Curvature 

\vx\  <  2299 

Cross-Range  Smear 

\vx\  <  25.53 

Range  Duplicates 

\vx\  <  0.7536 

Blind  Velocities 

vx  =  1V16.6  x  10-3 

Table  4:  Maximum  velocity  constraints  for  a  satellite  radar  system  with  the 
operating  parameters  given  in  Table  2.  For  the  blind  velocities,  N  denotes  any 
integer.  Velocities  are  normalized  by  the  platform  velocity. 


Term 

Limit 

Range  Walk 

\vy\  <  2.3 

Cross-Range  Duplicates 

\vy\  <  0.0435 

Range  Curvature 

Ur  <  106.2 

Cross-Range  Smear 

\vx\  <  2.0 

Range  Duplicates 

\vx\  <  0.1 

Blind  Velocities 

vx  =  N3.77  x  10~3 
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2.2  DPCA  matched  filter  mismatch 


The  previous  analysis  shows  the  output  of  the  matched  filter  when  the  filter  is  per¬ 
fectly  matched  to  the  incoming  data.  In  this  section  the  output  of  the  matched  filter 
when  there  is  an  error  in  the  parameters  (x0,  yo,  vx,  vy),  is  explored.  First  we  examine 
the  situation  when  the  matched  filter  is  misaligned  by  Xo  +  Ax.  The  resulting  output 
is  then  given  by: 


S(k,u,x)M(k,u,x  +  Ax)  =  exp  yj2(k  +  k0)  ^Ax  + 


exp 


(  . ,  ,  .  / .  dAx  Ax  9 

j2(k  +  k0)  \  Ax  + -^-u 


AT2  2A2 


u 


/  /  dAx  vyd  vxd  9  Ax  9N 

-  exp  \  j2(k  +  k0)  I  Ax  +  vxd  +  -^-u  -  —u  -  u  -  -^u 


dAx  v„d  vxd  2  Ax 


exp  j2(k  +  k0)  Ax  -  vxd  +  — —u  +  -A-u  +  -A- u 


A2  Xr  2A2  2A2 


u 


(  /,  ,  \  ( A  dAx  Ax  o\  \ 

=  2  exp  \  j2(k  +  fc0)  (Ax  +  u  -  )  I 

/  . ,  ,  /  A  ,  dAx  vvd  vxd  9  Ax  9N 

-  exp  f  j2(A  +  fc0)  (Ax  +  uxd  +  -  —  u  -  —u  -  —u 

( ,  .  ( .  ,  dAx  vvd  vxd  9  Ax 

-  exp  I  ]2(k  +  ko)  I  Ax  -  vxd  +  -p-n  +  -A-u  +  ~  2W“" 


(43) 


If  we  assume  that  the  radar  is  operating  under  the  five  criteria  listed  in  section  2.1 
(see  (36)  -  (40))  then  the  phase  terms  containing  and  f^u2  are  insignificant  and 
the  following  simplification  is  reached: 


f  (  d  /\  x  /\  x  ^ 

S(k,u,x)M(k,  u,x  +  Ax)  «  2  exp  \  j2{k  +  k0)  Ax  +  -  "■2 


-u 


X  2A2 


u 


—  exp  j2(k  +  ko)  Ax  +  vxd  + 


dAx  Ax 


A2 


-u  — 


2A2 


u 


exp  j2(k  +  ko)  Ax  —  vxd  + 


dAx  Ax 


A2 

C 


■U 


2A2 


u 


(44) 


Examination  of  the  phase  elements  of  (44)  reveals  that  the  output  of  the  matched 
filtering  operation  is  shifted  in  the  range  direction  by  Ax,  which  is  desirable  for 
imaging  purposes;  however,  the  filter  mismatch  also  introduces  linear  and  quadratic 
terms,  with  respect  to  u,  in  the  phase.  The  phase  term  k^A^-u  produces  range  walk 
in  the  resulting  SAR  image.  To  limit  the  range  walk,  the  shift  should  be  less  than 
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one  range  cell  at  the  limit  of  the  aperture.  This  constraint  produces  the  following 
inequality: 


07  d  |  Ax | 

^ fcm.ax  Umax  ^ 


AT2 


,  a  ,  7 rA2 

I  Ax  I  < 


dk'm.n.r.'U 


max  max 


(45) 


The  other  linear  phase  term,  koj^-u,  produces  a  shift  in  the  cross-range  dimension 
of  the  SAR  image.  This  does  not  affect  the  image  formation  process  of  the  moving 
target. 

The  phase  term  —j2(k  +  ko)X^u2  introduces  range  curvature  and  a  cross-range  smear 
component  in  the  resulting  image.  Limiting  the  effect  of  the  range  curvature  to  one 
resolution  cell  at  the  limit  of  the  aperture  requires: 


2  kr 


|Ax|  2 


2X2 


<  27T 


I  Ax  I  < 


2ttX2 


(46) 


k  ii 2 
A' max  max 


Limiting  the  cross-range  smearing  in  the  resulting  imagery  requires: 


2k, 


|Ax| 


7 r 


-7/“ 

2X2  max 


|Ax|  < 


t rX2 


4^0  V1  max 


(47) 


Next  we  examine  the  situation  when  the  filter  is  misaligned  by  y$  +  Ay.  In  this  case 
the  resulting  output  is  given  by: 


S(k,  u,  x)M (. k ,  u,  y  +  Ay) 


2  exp  ( j2 (k  +  k0)(^-~  ^ u 


Xr 


Xr 


dAy  Ay 

~X7  ~  X 


—  exp  I  j2{k  +  ho)  I  vxd  4 — —A  — -p-u - A- u  —  X^u2 


vyd 


vrd 


Xr 


2  Xr 


(  -o n  i  i  \  I  .1  ,  dAy  Ay  ,  vvd  ,  Lxd 

exp  \j2(k  +  ko)  I  -vxd  +  — - —  u  +  — u  + 

(48) 


Again,  assuming  that  radar  is  operating  such  that  the  phase  terms  and  ||4-u2are 
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insignificant  produces  the  following: 

S(k,u,x)M(k,u,y  +  Ay)  sa  2 exp  ^j2(k  +  k0) 

-  exp  ^j2(k  +  fc0)  (vxd  +  (49) 

-  exp  ^j2(k  +  ko)  vxd  +  . 

Examination  of  the  phase  of  elements  of  (49)  reveals  the  result  of  the  mismatch  pro¬ 
duces  not  only  an  offset  in  the  cross-range  direction  but  also  a  range  walk  component. 
The  offset  in  the  cross-range  dimension  has  no  effect  on  the  imaging  results.  However, 
limiting  the  range  walk  component  to  one  range  cell  at  the  extreme  of  the  aperture 
requires: 


2k  u 


max  y  max 
■**-c 


<  27 r 


i  a  i  7rXc 

<  k  v~ 

Ojm.ax 


max  max 


(50) 


For  the  purpose  of  illustration  the  tolerances  are  calculated  for  the  position  mis¬ 
matches  (Ax,  Ay  )  for  the  airborne  and  satellite  radar  systems  analyzed  in  Section 
2.1.  Table  5  and  Table  6  list  the  position  mismatch  constraints  for  the  airborne  and 
satellite  radar  geometry  respectively.  In  both  these  geometries  the  tolerance  for  the 
range  walk  and  the  cross-range  smear  are  the  most  dominant.  For  both  these  radar 
systems  the  matched  filter  would  have  to  be  recalculated  every  few  resolution  cells. 

Table  5:  Position  tolerances  of  matched  filter  for  airborne  radar  with  system 

parameters  given  in  Table  1. 


Term 

Limit 

Range  Walk 

|Ax  <  58.7  km 

Range  Walk 

A y\  <  29.37  m 

Range  Curvature 

Ax  <  2299  m 

Cross-Range  Smear 

Ax  <  6.38  m 

We  examine  the  situation  when  the  filter  is  mismatched  by  vx  +  Avx.  In  this  case  the 
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Table  6:  Position  tolerances  of  matched  filter  for  satellite  radar  with  system 


parameters  given  in  Table  2. 


Term 

Limit 

Range  Walk 

Ax  <  1959  km 

Range  Walk 

A y\  <  17.28  m 

Range  Curvature 

|Ax  <  796  m 

Cross-Range  Smear 

Ax  <  3.75  m 

S(k ,  u,  vx)M(k,  u,  vx  +  Avx)  =  exp  ^  j2(k  +  k0)  ^—dAvx  + 


+  exp 
—  exp 


dAvx  9  Avx  A 

“  'W“, 


vxd  2 


(j^k  +  k0)  [avxu  + 

^ j2{k  +  kQ)  (vxd  +  Avxu  -  ^ 

A  „ .  \  \ 


dAvx 


Vyd  UL-\UX 

u  +  ~xr« 


~2X? 


Avx  3 

“  -2Xf“ 


-dAv. 


exp  |  j2(k  +  h 

'U  +2X^U  ~  2Xf 


'x  ux 


1  A  Vvd 

vxd  +  Avxu  +  —i 

■A-C 


3A vxd  2  vxd  2  Au*,  3 
■  u 


2X? 


Again,  assuming  that  the  radar  is  operating  such  that  the  phase  terms  u 

-A.  c 

insignificant  produces  the  following: 


mi,  assuming  tiiat  tlie  r 
insignificant  produces 

S(k,  u,  vx)M{k ,  u,  vx  +  Avx)  «  exp  ^j2(k  +  k0)  ^-dAvx  +  A vxu  + 

^'2(/e  +  fc0)  ^Ar^w  + 

|j2(fc  +  fc0)  +  A vxu  + 

,  7.  dAvx  —  vxd  +  At>xn  -I 


am 


(j2(fc+ 

-  exp 

—  exp 

-  exp  j2(k  +  h 


3dAvx  9 
u 

c 

dA?;x  2  Aux  O 

— “  -  2X|“ 

dAvx  2 
■u  —  - 


Avxd  2  Auj;  3n 

-u 


2X1  U  2Xf 
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As  in  the  case  for  the  position  mismatch  the  velocity  mismatch  results  in  linear  and 
quadratic  phase  errors;  however  an  additional  third  order  phase  error  emerges.  The 
linear  and  quadratic  phase  errors  will  have  the  same  effects  as  discussed  previously. 
The  third  order  phase  error  term  introduces  higher  order  range  migration  (third 
order)  and  produces  asymmetrical  sidelobes  in  the  cross-range  dimension  of  the  SAR 
image  [8]. 

Limiting  the  effect  of  duplicate  targets  in  the  range  dimension  to  be  less  than  one 
range  cell  apart  requires: 


2‘kmaxd  27T 

(Ausl  < 


(53) 


h  rl 

rXjm.nx 


Limiting  the  range  walk  to  be  less  than  one  range  cell  at  the  extreme  limit  of  the 
aperture  requires  the  following  condition: 


2kmax  A'l’x  Umax  <  2tT 

|Anx|  < 


7 r 


(54) 


k max ^ 


max  max 


Limiting  the  range  curvature  to  be  less  than  one  range  cell  at  the  extreme  limit  of 
the  aperture  requires  the  following  condition: 


2  kri 


\Avx\d  9 

- 7 /  <7  27 T 

2,X^  Uj'max  ^ 


\Avx\  < 


2nX: 


(55) 


rlh  ii 2 
Ujn'maxu'max 


In  order  to  limit  the  cross-range  smearing  the  maximum  phase  variation  should  be 
less  than  7t/4  which  requires: 


<  - 


2A? 


I  An*  |  < 


ttAL2 


(56) 


4  dk0u 


2 

max 


Imposing  the  same  limit  on  the  third  order  range  migration  to  be  less  than  one  range 
resolution  cell  at  the  ends  of  the  aperture  produces  the  following  condition: 


2  krn 


I  An 


2A? 


x\  3 

tC 


<  27T 


|Adx|  < 


2nX: 


(57) 


h  n 3 

a'max  amax 
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Limiting  the  phase  error  to  be  less  than  7t/4  for  the  third  order  cross-range  smearing 
term  produces  the  following  condition: 


2  h 


\Avx 


3  n 

-V  <  — 

2X2  max  4 


|Avx|  < 


kouL 


(58) 


Table  7  and  Table  8  list  the  constraints  on  Avx  for  the  airborne  geometry  and  the 
satellite  geometry  respectively.  It  is  important  to  note  that  these  velocities  are  given 
relative  to  the  radar  platform  velocity.  In  both  these  geometries  the  limiting  factor 
is  driven  by  the  range  walk  and  third  order  cross-range  smearing  criterion.  This 
increased  sensitivity  will  result  in  a  large  number  of  matched  Liters  for  different  vx 
values  for  typical  radar  platform  velocities  and  expected  target  motion  velocities. 

We  examine  the  situation  when  the  filter  is  mismatched  by  vy  +  Avy.  In  this  case  the 


Table  7:  Avx  velocity  tolerances  for  matched  filter  for  airborne  radar  with  system 

parameters  given  in  Table  1. 


Term 

Limit 

Range  Duplicate 

\Avx\  <  0.75 

Range  Walk 

|Avx|  <  1.4  x  10"a 

Range  Curvature 

An^l  <  2298 

Cross-Range  Smear 

An^l  <  3261 

Third  Order  Range  Migration 

\Avx\  <  4.49 

Third  Order  Cross-Range  Smear 

IAuj.  <  49.9  x  10~3 

Table  8:  Avx  velocity  tolerances  for  matched  filter  for  satellite  radar  with  system 

parameters  given  in  Table  2. 


Term 

Limit 

Range  Duplicate 

|AUr|  <  0.1 

Range  Walk 

Avx  <  2  x  10~5 

Range  Curvature 

At^l  <  106.2 

Cross-Range  Smear 

Avx|  <  18472 

Third  Order  Range  Migration 

At^l  <  0.021 

Third  Order  Cross-Range  Smear 

\Avx\  <  4.0  x  10”4 
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output  of  the  matched  filter  is  given  by: 


S(k,  u,  vy)M (, k ,  u,  vy  +  A vy)  =  exp  j2(k  +  k0) 


[  2dAv„ 


An 


-u  — 


X, 


yu 2 


exp  j2(k  +  k0) 


/  ClAVy 


An 


yu 2 


-  exp  j2(k  +  ko)  vxd  + 


exp  j2(k  +  fc0)  ~vxd  + 


Xc  ,  , 

dAVy  Vyd  A  Vy  9  Vj-d. 

~x7u~ircu~ic;u  ~2xzu\ 

2dAv 


Vyd  A  Vy  9  VXd  9  ' 

Y~u  +  lc.u  ~  ~x7u  +  ~2JQU  , 


(59) 


Again,  assuming  that  the  radar  is  operating  under  the  conditions  given  by  (36)  - 
(40)  then  the  phase  terms  ^ u  and  f^u2  are  insignificant  produces  the  following 
simplification: 


S(k,  u,  vy)M(k ,  u,  vy  +  Avy)  =  exp  \  j2(k  +  k0) 

+  exp  j2(k  +  k0 


f  2dAv. 


{  xcV 

(  dAv. 


Av. 


u 


Xc 

An 


yu 2 


V  Ac 


u  — 


X, 


yu 2 


exp  j2(k  +  k0)  vxd  + 


dAv,, 


Av . 


Xr 


-u 


Xr 


yu2 


exp  j2(k  +  k0)  -vxd  + 


2dAvy 


Av. 


Xr 


-u 


Xr. 


yu 2 


(60) 


Examination  of  the  phase  elements  of  (60)  reveals  the  presence  of  undesirable  linear 
and  quadratic  phase  terms.  The  phase  terms  produce  effects  similar  to  those  previ¬ 
ously  analyzed.  Limiting  the  range  walk  to  be  less  than  one  range  cell  at  the  extreme 
limit  of  the  aperture  produces  following  condition: 


2k  r 


| Avy\  d 


Xr 


Ur, 


<  2n 


\AVy\  K 


7 rXr 


(61) 


dkmaxumax 


Limiting  the  range  curvature  to  be  less  than  one  range  cell  at  the  extreme  limit  of 
the  aperture,  results  in  the  following  condition: 


2kr. 


I  An 


A? 


y'u 2 


<  2n 


I  A  I  nXC 

|Auy|  < 


(62) 


k  ?/ 2 

r^max  umax 
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Limiting  the  phase  error  to  be  less  than  7t/4,  in  order  to  limit  the  cross-range  smearing, 
requires  the  following  condition: 


2kJXAu* 


X2 

c 


7 r 

max  ^  ^ 


|Auy|  < 


7rX2 


8k0u2n 


(63) 


Table  9  and  Table  10  list  the  constraints  on  Avy  for  the  airborne  and  satellite  radar 
geometries  respectively  .  ft  is  important  to  note  that  these  velocities  are  given  rel¬ 
ative  to  the  radar  platform  velocity.  For  both  these  geometries  the  most  restrictive 
criterion  is  related  to  the  range  walk.  However,  for  the  typical  radar  platform  veloci¬ 
ties  and  expected  range  of  target  motion  velocities  these  restrictions  are  not  expected 
to  be  prohibitive.  From  the  analysis  of  the  tolerances  for  both  position  and  velocity 
mismatches  in  the  matched  filter  it  is  apparent  that  the  dominating  factors  are  the 
mismatches  in  both  range  and  cross-range  positions  and  the  target’s  velocity  in  the 
range  direction,  i.e.,  vx. 

Table  9:  Avy  velocity  tolerances  for  matched  filter  for  airborne  radar  with  system 

parameters  given  in  Table  1. 


Term 

Limit 

Range  Walk 

\AVy 

<  29.3 

Range  Curvature 

\AVy 

;  1149 

Cross-Range  Smear 

l  AVy 

<  1630 

Table  10:  Avy  velocity  tolerances  for  matched  hlter  for  satellite  radar  with  system 

parameters  given  in  Table  2. 


Term 

Limit 

Range  Walk 

|Auy  <  2.3 

Range  Curvature 

\AVy 

|  <  398 

Cross-Range  Smear 

|AuJ 

<  69272 

22 


DRDC-RDDC-2014-R51 


3  Summary  and  conclusions 


This  report  develops  the  matched  filter  necessary  for  focusing  a  moving  target  signa¬ 
ture  after  DPCA  processing.  Due  to  the  signal  subtraction  in  the  DPCA  processing, 
the  resulting  imagery  will  contain  artifacts.  The  appropriate  limits  on  the  range  and 
cross-range  velocities  are  derived  to  limit  these  artifacts  and  the  corresponding  values 
calculated  for  typical  airborne  and  satellite  geometries.  In  both  these  geometries  the 
most  stringent  constraint  is  imposed  by  the  limitation  on  the  range  and  cross-range 
duplicates,  which  limit  the  range  and  cross-range  velocity  respectively  of  the  target 
of  interest. 

Additionally,  the  effect  of  filter  mismatch  to  the  moving  target’s  position  and  velocity 
has  also  been  analyzed  and  the  corresponding  limits  on  the  targets  position  and 
velocities  used  in  the  matched  filter  have  been  derived.  For  both  the  airborne  and 
satellite  geometries  the  limitation  on  cross-range  smearing  in  the  imagery  imposes 
the  tightest  constraint  on  the  range  mismatch  to  the  target.  The  mismatch  in  the 
target’s  range  velocity  is  most  tightly  constrained  by  the  range  walk  criterion  for  the 
both  airborne  and  satellite  geometries.  However,  the  effect  on  third  order  cross-range 
smearing  should  also  be  examined.  For  both  the  airborne  and  satellite  geometries 
the  mismatch  in  the  target’s  cross-range  velocity  is  most  tightly  constrained  by  the 
limitation  on  range  walk.  However  due  to  the  relative  insensitivity  of  the  image 
formation  process  to  cross-range  velocity  mismatch,  this  limitation  is  unlikely  to  pose 
any  problems  for  most  SAR  systems. 
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Annex  A:  Taylor  series  approximations  for 

cos  02  and  sin  02 


In  this  section  we  derive  the  Taylor  series  formulas  for  cos  02  and  sin  02.  These  are 
the  most  general  cases.  Other  simpler  cases  can  be  found  by  setting  the  appropriate 
values  to  zero  e.g.,  Yc  =  0.  Further  approximations  can  also  be  made  by  assuming 
d  =  0,  or  d  -C  Xc. 

A.1  Approximation  for  cos  02{u) 

From  (26)  the  definition  of  cos  62  (u)  is  given  as: 


cos  62(11)  =  .  (A.l) 

y/Xl  +  iY'-U  +  d)* 

Evaluating  cos  6*2  (u)  at  u  =  0  gives: 

cos  #2(0)  =  7  (A. 2) 

v'.V2  +  (V  +  df 

Taking  the  first  derivative  of  cos 02{u)  gives: 


d  cos  62  (u) 
du 


-\xc  (xc2  +  (Yc  -  u  +  df)  1  (—2  ){YC  - u  +  d ) 


XC(YC  -u  +  d) 


(A.3) 


(Xc2  +  (Yc-u  +  d)2)* 

Evaluating  the  first  derivative  at  u  =  0  gives: 


d  cos  0o(u) 


du 


XC(YC  +  d) 


,=o  {X2  +  {Yc  +  d)2y 


(A.4) 


Taking  the  second  derivative  of  cos02(u)  gives: 


d 2  cos  62  (u) 
du2 


_ Xc _ 

{X2  +  {Yc-u  +  d)2)i 
_ Xc _ 

{X2  +  (Yc-u  +  d)2)^ 


3  XC(YC  —  u  +  d) 

2  (Xc2  +  (Yc-u  +  d)2f* 
Xc{Yc-u  +  d )2 
(Xc2  +  (Yc-u  +  d)2)^ 


(-2  ){Yc-u  +  d) 


(A.5) 
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Evaluating  the  second  derivative  at  u  =  0  gives: 


d 2  cos 9o(u) 


du 2 


Xr. 


+  3- 


XC(YC  +  df 


i=0  (Xc2  +  (Yc  +  d)2)*  (Ac2  +  (yc  +  d)2)2 


(A.6) 


The  Taylor  series  approximation  for  cos  92  (u)  is  then  given  as: 

AT  AC(YC  +  d) 


COS  6*2  (u) 


+ 


\jxc2  +  (Yc  +  d)2  (A c2  +  (yc  +  d)2)3 


Ac(yc  +  d)s 


AT 


(AT2  +  (yc  +  d)2)5  (ac2  +  (yc  +  d)2y 

For  the  broadside  case,  set  Yc  =  0  which  gives: 


u 

Y 


COS  6*2  (u)  = 


AT 


+ 


AT  d 


;U  +  I  3- 


X rd2 


AT 


VAC2  +  d2  (x2  +  d2Y  \  (Ac2  +  d2)§  (Ac2  +  d2)§ 
d  (  d2  1  \  u2 

+  x2U+{3x?~x2)y- 

In  the  case  where  d  -C  AT,  and  0  -C  Xc  then: 

„  ,  .  d  u2 

cos 92{u)  «  1  +  —  u  - 


AT2  2  A2  ‘ 


In  the  particular  case  when  d  =  0  then: 

cosd(w)  ~  1 


u 


2A2 


(A.7) 


u 


(A.8) 


(A.9) 


(A.10) 


A.2  Approximation  for  sin  02 (w) 

From  (27)  the  definition  of  sin  #2  (it)  is  given  as: 

sin02(w)  =  6  .  (A. 11) 

\J  X2  +  (Yc  —  u  +  d)2 

Evaluating  sin02(w)  at  u  =  0  gives: 

sin#2(0)  =  ,  +  ^  (A. 12) 

\jX2  +  (Yc  +  d)2 
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Taking  the  first  derivative  of  sin  62  (u)  gives: 
d  sin  d2  (u)  1  1  (Yc  —  u  +  d) 


du 


^X2  +  (Yc-u  +  d)2  2  (xc2  +  (Yc  -  u  +  d)rf 

1  (Yc  —  u  +  d)2 


(-2  ){Yc-u  +  d) 


{X2  +  (Yc  -  u  +  d)2Y  +  (Xc2  +  (Yc  -  u  +  d)2)^ 


Evaluating  the  first  derivative  at  u  =  0  gives: 


d  sin  92  (u) 


du 


+ 


(Yc  +  d)‘ 


i—0  (Xc*  +  (Y)  +  <i)2)’  (AV  +  (Y)  +  d)?)-’ 


(A-13) 


(A. 14) 


The  Taylor  series  approximation  for  sin  62  (u)  is  then  given  as: 


sin  62  (u) 


Yc  +  d 


+ 


(Yc  +  d)s 


^X2  +  (Yc  +  d)2  \(Xc2  +  (Yc  +  d)2f*  (Xc2  +  (Yc  +  d)2)^ 


u. 


(A. 15) 


For  the  broadside  case,  set  Yc  =  0  which  gives: 


sin  6*2  (u)  = 


d 


+ 


d2 


VXj+d2  \(Xc2  +  d2)5 

In  the  case  where  d  -C  Xc,  and  0  -C  Xc  then: 

•  a  (  \  d  u 
Sill  t) 2{U)  ~  ~zz - — ■ 

In  the  particular  case  when  d  —  0  then: 


(x2  +  d2y\ 


u. 


(A. 16) 


AT  Ac 


sin0(u)  =  — 


u 

AT 


(A. 17) 


(A. 18) 
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